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Abstract 

We  present  experimental  flux-temperature  phase  diagrams  for  surface  reconstruction  transitions  on  the  6.1  A  com¬ 
pound  semiconductors.  The  phase  transitions  occur  within  or  near  typical  substrate  temperature  ranges  for  growth  of 
these  materials  by  molecular  beam  epitaxy  and  therefore  provide  a  convenient  temperature  standard  for  optimizing 
growth  conditions.  Phase  boundaries  for  InAs  (00  1)  [(2  x  4)  ->  (4  x  2)],  AlSb  (00  1)  [c(4  x  4)  -» (1  x  3)],  and  GaSb  (00  1) 
[(2  x  5)  ->  (1  x  3)]  are  presented  as  a  function  of  substrate  temperature  and  Group  V-limited  growth  rate  (proportional  to 
flux),  for  both  cracked  and  uncracked  Group  V  species.  We  discuss  differences  between  materials  in  the  slopes  and  offsets 
of  the  phase  boundaries  for  both  types  of  Group  V  species.  ©  2000  Elsevier  Science  B.V.  All  rights  reserved. 

PACS:  18.15.Hi;  81.05.Ea;  68.35.Bs;  68.35.Rh 

Keywords:  Molecular  beam  epitaxy;  Gallium  antimonide;  Aluminum  antimonide;  Indium  arsenide;  Surface  reconstruction 


1.  Introduction 

The  compound  semiconductors  with  6.1  A  lattice 
constant,  including  InAs,  GaSb,  and  AlSb,  have 
received  increasing  attention  in  recent  years  be¬ 
cause  of  their  potential  for  electronic  and  optical 
device  applications  [1-6].  A  unique  combination  of 
band  offset  (broken  gap  type  II,  between  InAs  and 
GaSb),  high  InAs  electron  mobility,  and  InAs  sur¬ 
face  Fermi  level  pinning  permit  novel  device 
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designs  and  improved  performance  that  are  diffi¬ 
cult  to  achieve  with  other  compound  semiconduc¬ 
tors.  In  order  to  realize  these  benefits,  epitaxial 
growth  and  device  processing  technologies  for  the 
6.1  A  compound  semiconductors  must  advance  to 
the  level  of  more  mature  material  systems  such  as 
GaAs/AlGaAs  and  InP/InGaAs. 

In  order  to  grow  the  highest  quality  compound 
semiconductor  thin  Aims  by  molecular  beam  epi¬ 
taxy  (MBE),  precursor  fluxes  and  substrate  temper¬ 
ature  should  be  optimized  [7,8].  These  growth 
parameters  affect  the  rates  of  adsorption,  desorp¬ 
tion,  and  surface  diffusion  of  Group  III  and  V 
species,  and  therefore  influence  the  quality  of  crys¬ 
tal  growth.  Activation  barriers  and  preexponentials 
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for  these  processes  depend  on  the  structure  and 
stoichiometry  of  the  material  surface  (see  e.g.  Refs. 
[9,10]),  which  are  in  turn  affected  by  temperature 
and  beam  flux.  To  achieve  the  best  growth  condi¬ 
tions,  it  is  therefore  desirable  to  know  the  recon¬ 
struction  phase  boundaries  over  typical  beam  flux 
and  substrate  temperature  ranges. 

The  flux-temperature  phase  diagram  is  also 
convenient  for  estimating  substrate  temperature, 
particularly,  when  standard  techniques  for 
temperature  measurement  are  unavailable.  Remote 
thermocouple  measurements  in  vacuum  are  notori¬ 
ously  inaccurate.  Furthermore,  their  calibration 
can  change  during  growth  of  small-bandgap 
epilayers  (e.g.  InAs  and  GaSb),  due  to  radiative 
heating  that  increases  with  film  thickness  [11,12]. 
The  accuracy  of  optical  pyrometry  is  limited  by 
stray  light  from  the  substrate  heater  or  source 
ovens,  coating  of  the  viewport,  changes  in  substrate 
emissivity,  and  limited  sensitivity  below  500°C 
[13].  Transmission  thermometry,  which  monitors 
the  shift  in  the  optical  transmission  edge  of  the 
substrate,  is  an  effective  alternative  to  pyrometry 
[14-16],  However,  this  technique  is  not  possible  for 
InAs  substrates  due  to  free  electron  absorption,  and 
measurements  for  6.1  A  heteroepitaxy  on  GaAs 
substrates  are  not  possible  when  a  small-bandgap 
epilayer  is  thick  enough  to  absorb  most  of 
the  source  light  ( >  1  pm).  When  these  or 
other  techniques  are  not  available,  observation 
of  known  surface  reconstruction  transitions 
remains  a  convenient  means  to  estimate  substrate 
temperatures. 

There  have  been  several  reports  of  reconstruc¬ 
tion  phase  diagrams  for  III-V  materials,  especially 
for  GaAs.  Early  work  on  reconstruction  transitions 
[17,18]  includes  measurements  for  all  zincblende 
Ill-V  compounds  at  a  fixed  Group  V  flux  [18]. 
Full  phase  diagrams  have  been  reported  for  GaAs 
(001)  [19-23],  GaAs  (111)  [24],  AlAs  (001) 
[19,20,25],  AlGaAs  (001)  [22],  and  for  the  6.1  A 
materials  [16,23,26,27],  although  not  all  with  calib¬ 
rated  fluxes  of  both  cracked  and  uncracked  Group 
V  species.  Recent  results  on  GaSb  from  our  group 
were  presented  previously  [16]  and  are  discussed  in 
this  work  as  well. 

Phase  diagrams  typically  describe  thermodyn¬ 
amic  equilibrium  conditions.  However,  for  MBE 


conditions  during  surface  stabilization  by  a  Group 
V  flux,  the  system  is  not  in  strict  equilibrium:  the 
vapor  and  solid  phases  are  at  different  temper¬ 
atures;  the  vapor  velocity  distribution  is  aniso¬ 
tropic;  the  incident  and  desorbing  species  may  not 
be  the  same.  Nevertheless,  equilibrium  thermo¬ 
dynamics  may  be  used  as  an  approximate  descrip¬ 
tion  of  a  Group  V-stabilized  surface  [7,8].  We 
adopt  this  view,  although  in  Section  4  we  discuss 
some  evidence  for  non-equilibrium  behavior. 

2.  Experimental  method 

The  6.1  A  semiconductors  studied  in  this  report 
were  grown  by  solid-source  MBE.  Growth  rates 
and  surface  reconstruction  symmetry  were  moni¬ 
tored  with  reflection  high-energy  electron  diffrac¬ 
tion  (RHEED).  A  10  kV  electron  beam  was  incident 
on  the  sample  surface  between  1.5°  and  2.5°.  The 
beam  had  been  previously  adjusted  to  focus  on  the 
RHEED  phosphor  screen  (~70cm  from  electron 
gun  to  screen)  without  the  sample  in  place.  All  three 
6.1  A  materials  were  grown  on  semi-insulating 
GaAs  (0  0 1)  substrates,  and  the  GaSb  measurement 
was  also  repeated  on  GaSb  substrates,  producing 
no  shift  of  the  phase  boundary  [16].  AlSb  and 
GaSb  layers  on  GaAs  were  roughly  1  pm  thick, 
while  InAs  was  grown  on  a  1  pm  AlSb  buffer.  Thin 
layers  of  InAs  (1000  A)  were  used  in  order  to  min¬ 
imize  absorption  of  radiation  from  the  substrate 
heater,  which  is  needed  for  transmission  thermom¬ 
etry.  In  order  to  confirm  that  strain  in  the  thin  InAs 
layer  (1.3%  lattice  mismatch  with  AlSb)  does  not 
significantly  influence  the  reconstruction  transition 
temperature,  we  measured  this  temperature 
for  InAs  (001)  [(2  x  4) ->  (4  x  2)]  under  oa 

background  arsenic  flux  for  thickness  up  to  8000  A. 
No  dependence  on  the  InAs  layer  thickness  was 
observed. 

Substrate  temperatures  were  measured  by  trans¬ 
mission  thermometry.  This  technique  relies  on  the 
change  in  the  substrate  bandgap  with  temperature 
[14-16],  Briefly,  light  from  the  substrate  heater  is 
partially  transmitted  through  the  substrate  and 
exits  the  chamber  through  a  viewport.  The  light  is 
chopped  (for  phase-sensitive  detection),  focused 
into  an  optical  fiber,  dispersed  by  a  grating 
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spectrometer,  and  finally  detected  with  a  photo¬ 
diode.  The  transmission  spectrum  of  the  substrate 
is  obtained  by  taking  the  ratio  of  the  spectrum 
obtained  with  the  substrate  in  place  to  a  back¬ 
ground  spectrum  of  the  bare  substrate  heater  emis¬ 
sion.  For  GaSb  substrates,  the  relation  between 
transmission  edge  and  substrate  temperature  was 
taken  from  Ref.  [16].  In  the  appendix,  we  present 
a  new  calibration  for  GaAs  substrates,  in  order  to 
update  deficiencies  in  an  earlier  version. 

EPI  cracker  models  500V-As  and  175-Sb  were 
used  to  produce  cracked  Group  V  species.  By 
choosing  the  appropriate  operating  conditions  for 
the  arsenic  source,  it  is  straightforward  to  produce 
either  arsenic  dimers  (As2)  or  tetramers  (As4),  but 
the  antimony  cracker  source  (cracker  temperature 
900°C)  produces  a  mixture  of  Sb2  and  Sb  atoms 
due  to  the  small  Sb2  bond  strength.  A  detailed 
characterization  (using  the  175-Sb)  was  reported  in 
Ref.  [28]  and  showed  that  the  source  produces 
a  large  flux-dependent  proportion  of  Sb  atoms  un¬ 
der  operating  conditions  typical  for  our  work.  Anti¬ 
mony  dimer  beams  discussed  in  later  sections  of 
this  paper  should  therefore  be  understood  to  in¬ 
clude  atoms  as  well.  We  use  a  separate  antimony 
source  (operating  temperatures  <500°C)  to  pro¬ 
duce  tetramers. 

Although  molecular  beam  fluxes  are  commonly 
reported  as  ion  gauge  pressure  readings,  such  read¬ 
ings  are  often  not  reproducible  between  laborator¬ 
ies  due  to  differences  in  ion  gauge  geometry, 
electrode  voltages,  filament  and  electrode  coatings, 
etc.  Furthermore,  the  conversion  from  an  ion  gauge 
reading  to  beam  fluxes  requires  quantitative  values 
for  the  species-dependent  ionization  cross-sections. 
Because  of  these  shortcomings,  we  also  present 
phase  diagrams  with  the  flux  expressed  in  units  of 
growth  rate  at  a  11I/V  flux  ratio  for  which  the 
Group  V  flux  is  limiting.  This  approach  has  been 
used  previously  in  studies  of  epitaxial  growth  and 
surface  reconstructions  on  GaAs  (see  e.g.  Refs. 
[16,29]). 

We  used  the  following  procedure  to  convert  to 
Group  V-limited  growth  rate  units.  These  measure¬ 
ments  are  made  separately  from  the  phase 
transition  measurements.  For  a  fixed  Group  V  flux, 
the  growth  rate  is  measured  from  RHEED  oscilla¬ 
tions  as  a  function  of  Group  111  cell  temperatures. 


835  840  845  850 

indium  cell  temperature  (°C) 

Fig.  1.  Typical  Group  V-limited  growth  rate  measurement,  here 
for  InAs  with  As2  (3.2  x  10~6Torr  beam  equivalent  pressure). 
Arsenic  flux  limits  growth  at  0.802  ML/s,  yielding 
g2  =  4.0  x  10“  6  Torr/(ML/s). 


As  the  Group  111  flux  is  increased,  the  growth  rate 
eventually  saturates  as  growth  becomes  limited  by 
the  Group  V  flux.  An  example  of  this  behavior  is 
shown  in  Fig.  1  for  InAs.  Dividing  the  Group  V  flux 
(beam-equivalent  pressure  units,  torr)  by  the  satu¬ 
ration  growth  rate  (ML/s)  gives  the  values  g2  (for 
dimers)  and  g4  (for  tetramers)  that  are  used  to 
rescale  the  vertical  axes  of  the  phase  diagrams. 
A  Group  Ill-rich  reconstruction  RHEED  pattern  is 
often  observed  during  the  V-limited  growth  oscilla¬ 
tions.  We  found  that  g2  and  g4  were  independent  of 
Group  V  flux,  with  the  only  exception  at  the  lowest 
fluxes  of  arsenic,  where  non-beam  arsenic  species 
may  contribute  a  significant  fraction  to  the  ion 
gauge  reading.  The  calibration  was  carried  out  at 
420°C,  430°C,  and  520°C  for  InAs,  GaSb,  and  AlSb, 
respectively.  These  temperatures  should  be  roughly 
reproduced  when  repeating  the  calibration,  because 
the  Group  V  incorporation  coefficients  during 
growth  on  (00  1)  surfaces  begin  to  decrease  at  high¬ 
er  substrate  temperature  (see  e.g.  Refs.  [29,3c]1  for 
typical  behavior). 


1  The  InAs  phase  diagram  in  Ref.  [23]  for  the  (2  x  4)  ->  (4  x  2) 
transition  during  growth  should  correspond  roughly  to  the 
temperature  dependence  of  V-limited  growth  as  defined  by  gx  in 
this  work. 


A  S.  Bracker  et  al.  /  Journal  of  Crystal  Growth  220  (2000)  384-392 


387 


In  As 

GaSb 

AlSb 

[110] 

[110] 

[110] 

(4x2) 

(1x3) 
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Fig.  2.  Schematic  representation  of  RHEED  pattern  evolution 
as  substrate  temperature  is  ramped  through  phase  transitions. 
Stars  mark  ihe  points  chosen  as  empirical  definitions  of  the 
phase  transitions. 


3.  Phase  diagrams 

Our  experimental  definitions  for  the  reconstruc¬ 
tion  phase  transitions  are  based  on  changes  in  the 
RHEED  pattern  (Fig.  2).  Our  primary  intention  in 
choosing  these  definitions  was  that  the  measure¬ 
ments  could  be  reproduced  in  other  laboratories. 
We  are  not  concerned  with  whether  the  corre¬ 
sponding  changes  in  surface  structure  satisfy  more 
rigorous  criteria  for  a  phase  transition.  Changes  in 
surface  periodicity  occur  at  different  temperatures 
for  different  crystal  directions,  and  mixed  phases  or 
domains  may  exist  around  the  transition  temper¬ 
ature.  We  simply  represent  the  transition  region  as 
a  line,  while  recognizing  that  additional  surface 
structures  may  exist. 

Several  sources  of  experimental  error  contribute 
to  uncertainty  in  the  phase  transition  temperatures. 
The  two  most  important  errors  in  our  experience 
were  the  position  of  the  RHEED  electron  beam  on 
the  substrate  (due  to  lateral  temperature  gradients) 
and  unsteady  ramp  rates.  A  background  pressure  of 
undesired  Group  V  species  is  likely  to  shift  the 
transition  temperature,  especially  for  antimonides 
when  arsenic  is  a  background  species.  We  kept  the 
arsenic  source  at  idle  temperature  during  antimon- 
ide  measurements,  and  vice  versa,  but  unfortunate¬ 
ly  this  may  not  be  possible  during  heterostructure 
growths  that  require  both  Group  V  species.  Finally, 
there  is  some  subjectivity  in  judging  changes  in 
diffraction  streaks  as  described  in  Fig.  2.  Given 


Q. 


450  440  430  420  410  400  390  380 


Substrate  Temperature  (°C) 

Fig.  3.  Reconstruction  phase  diagrams,  with  Group  V  flux  ex¬ 
pressed  in  beam-equivalent  pressure  units  (s2  and  s4  in  Torr). 
Horizontal  axis  is  linear  in  reciprocal  absolute  temperature. 
Solid  lines  are  fits  using  the  expression  and  parameters  in 
Table  1. 


these  errors,  repeated  measurements  rarely  varied 
by  more  than  10°C. 

Figs.  3  and  4  show  the  phase  diagrams  for  surface 
reconstructions  on  In  As,  GaSb,  and  AlSb.  We  used 
a  simple  exponential  form  to  fit  the  data  points,  as 
described  in  Table  1.  The  vertical  axes  of  Fig.  3  are 
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Substrate  Temperature  (°C) 


Fig.  4.  Reconstruction  phase  diagrams,  with  Group  V  flux  ex¬ 
pressed  in  units  of  Group  V-limited  growth  rate  (s2/gi  and  s4/g4 
in  ML/s).  Ftorizontal  axis  is  linear  in  reciprocal  absolute  temper¬ 
ature.  Solid  lines  are  fits  using  the  expression  and  parameters  in 
Table  1. 


expressed  in  the  units  of  beam-equivalent  pressure 
(Torr)  as  measured  with  an  ion  gauge,  which  is 
proportional  to  the  absolute  Group  V  atomic  flux 
required  to  stabilize  the  surface  at  the  phase  bound¬ 
ary.  The  plots  on  this  scale  are  referred  to  sub¬ 
sequently  as  s2  and  s4,  for  stabilization  under 


Table  1 

Parameters  used  in  the  expression  flux  =  exp  (  —  E/kB  T  +  B) 
for  fits  in  Figs.  3  and  4.  BTorr  and  BML/S  give  flux  in  units  of 
beam-equivalent  pressure  (Torr)  and  Group  V-limited  growth 
rate  (ML/s),  respectively 


E  (eV) 

^Torr 

^ML/s 

InAs 

As2 

4.16 

47.5 

60.0 

As4 

2.95 

31.5 

43.5 

AlSb 

Sb2 

4.64 

62.6 

76.4 

Sb4 

3.10 

39.0 

52.5 

GaSb 

Sb2 

1.78 

16.3 

30.1 

Sb4 

1.72 

15.7 

29.2 

a  Group  V  flux  of  dimers  or  tetramers,  respectively. 
This  form  of  the  phase  diagram  may  be  used  as 
a  quick  reference  for  estimating  substrate  temper¬ 
ature,  but  these  values  may  not  compare  well  be¬ 
tween  laboratories,  for  the  reasons  discussed 
earlier. 

Fig.  4  shows  the  phase  diagrams  with  vertical 
axes  rescaled  in  units  of  V-limited  growth  rate 
(s2/g2  and  s4/t/4).  The  divisors  g2  and  g4  (for 
growth  with  dimers  and  tetramers,  respectively)  are 
proportional  to  the  incident  Group  V  flux  (Torr) 
during  V-limited  growth  divided  by  the  Group 
V-limited  growth  rate  in  ML/s.  They  are 
t/2(AlSb)  =  1.0  x  10“6,  ry4( AlSb)  =  1.3  x  1CT6, 

t/2(GaSb)  =  1.0  xl(T6,  g4(GaSb)  =  1.3  x  10-6, 
t/2(InAs)  =  3.9  x  10-6,  g4(InAs)  =  6.4  x  10-6,  with 
units  of  Torr/(ML/s).  The  units  sx/gx  are  still  pro¬ 
portional  to  the  incident  Group  V  flux  required  for 
stabilization,  but  they  are  more  easily  transferred 
between  laboratories  than  sx,  because  ion  gauge 
sensitivities  are  divided  out.2  We  stress  that  gx  was 
measured  separately  and  at  only  one  substrate  tem¬ 
perature  for  purposes  of  rescaling  each  phase 


2  When  gx  and  sx  are  both  measured  in  the  temperature  range 
where  gx  is  independent  of  temperature,  then  the  units  sx/gx  are 
also  equal  to  the  rate  at  which  Group  V  atoms  would  be 
incorporated  into  the  crystal  during  a  V-limited  growth 
measurement  for  flux  sx.  This  condition  is  approximately  true 
for  GaSb  and  AlSb.  However  for  InAs,  the  reconstruction  phase 
boundaries  were  measured  at  temperatures  far  higher  than 
where  gx  was  measured  and  in  a  range  where  gx  rises  rapidly 
with  temperature. 
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boundary  in  Fig.  3  to  produce  Fig.  4.  No  growth 
occurred  during  the  phase  transition  measure¬ 
ments.  The  slopes  of  the  phase  boundaries  reflect 
the  thermodynamics  and  kinetics  of  stabilization 
only. 

The  changes  in  the  RHEED  pattern  used  to 
identify  the  phase  transitions  are  marked  by  stars  in 
Fig.  2.  For  all  measurements,  the  substrate  temper¬ 
ature  was  ramped  upward  at  a  rate  of  2°  per  min¬ 
ute.  For  AlSb  (001)  [c(4  x  4)  ->  (1  x  3)], 3  we  defined 
the  transition  temperature  by  equal  intensities  for 
the  1/3  and  2/3  streaks  of  the  3  x  pattern  and  the 
1/2  streak  of  the  2  x  pattern  [from  c(4  x  4)],  in  the 
[110]  direction. 

The  InAs  (001)  [(2  x  4)  ->  (4  x  2)]  transition  was 
defined  as  the  first  appearance  of  the  4  x  diffraction 
streaks  in  the  [110]  direction.  Yamaguchi  and 
Florikoshi  have  also  studied  this  transition  in  detail 
for  As4  flux  [26].  They  used  a  specialized  apparatus 
to  monitor  changes  in  the  RHEED  specular  inten¬ 
sity  and  measured  temperature  and  pressure  with 
pyrometry  and  an  ion  gauge,  respectively.  Their 
measured  phase  boundary  for  As2  flux  occurred 
roughly  30°C  lower  than  ours,  with  flux  expressed 
as  an  ion  gauge  reading.  This  discrepancy  is  prob¬ 
ably  a  consequence  of  a  different  experimental  def¬ 
inition  for  the  transition,  and  possibly  also  due  to 
inaccuracies  associated  with  pyrometer  or  ion 
gauge  measurements. 

The  GaSb  [(2  x  5)  ->  (1  x  3)]  phase  diagram  was 
recently  reported  by  our  group  [16].4  The  phase 
diagram  from  that  work  is  reproduced  here  in  Figs. 
3  and  4.  The  transition  was  defined  by  the  appear¬ 
ance  of  a  pseudo-(l  x  3)  RHEED  pattern,  as  the  1/5 
and  4/5  streaks  disappear,  and  before  the  2/5  and 
3/5  streaks  move  apart  into  the  typical  positions  for 
the  (1x3)  pattern. 

These  definitions,  combined  with  the  experi¬ 
mental  techniques  described  above  for  obtaining 


3  Recent  work  (Ref.  [40])  has  shown  that  the  AlSb  (1  x  3)  is  in 
fact  a  (4x3)  structure,  but  the  4  x  periodicity  is  usually  not 
observed  in  RHEED. 

4  The  GaSb  “(2  x  5)”  reconstruction  discussed  here  appears  in 
scanning  tunneling  microscope  images  as  (2  x  10)  or  c(2  x  10) 

unit  cells  (see  Ref.  [41]).  In  RHEED,  a  (1  x  5)  or  (2  x  5)  pattern  is 
observed. 


the  flux  scale  in  Fig.  4,  should  allow  reproducible 
measurement  of  the  phase  boundaries  in  other  la¬ 
boratories.  Because  the  phase  boundaries  fall  in  or 
near  the  usual  growth  temperature  ranges  for  the 
6.1  A  materials,  they  should  prove  useful  for  tem¬ 
perature  calibration  when  other  techniques  are  in¬ 
accurate  or  unavailable.  Although  the  AlSb  (001) 
[c(4  x  4)  — ►  ( 1  x  3)]  transition  falls  considerably  be¬ 
low  the  optimal  temperature  range  for  AlSb  growth 
(500-600°C),  most  device  structures  containing 
AlSb  also  contain  InAs  and  GaSb,  which  are  grown 
at  temperatures  below  500°C.  The  GaSb  (001) 
[(2  x  5)  ->  (1  x  3)]  transition  falls  within  the  com¬ 
mon  growth  temperature  range  for  InAs  and  GaSb 
(400-500°C),  while  the  InAs  (001)  [(2x4)-> 
(4  x  2)]  transition  lies  above  this  range. 

We  next  discuss  the  significant  differences  be¬ 
tween  materials  in  the  offsets  and  slopes  of  the 
phase  boundaries  for  cracked  and  uncracked 
Group  V  species. 

4.  Reactivity  of  Group  V  dimers  and  tetramers 

The  reactivity  of  uncracked  (tetramer)  and 
cracked  (dimer5)  Group  V  species  on  III-V  surfaces 
has  long  been  an  area  of  fundamental  and  practical 
interest.  The  form  of  the  incoming  species  will  affect 
the  rates  of  adsorption,  and  quite  possibly  diffusion 
and  desorption  as  well.  Adsorption  of  gas  phase 
Group  V  dimers  may  be  as  simple  as  a  bond  length 
change  and  rearrangement  on  the  surface  to  form 
a  surface  dimer.  In  contrast,  the  tetrahedral  Group 
V  tetramers  must  undergo  a  more  complicated 
electronic  and  nuclear  rearrangement  to  produce 
surface  dimers,  resulting  in  slower  rates.  Some 
experimental  results  for  gallium  arsenide  growth 
suggest  that  As4  chemisorption  requires  the 
interaction  of  two  neighboring  physisorbed 


5  For  cracked  antimony,  this  discussion  assumes  that  the 
monomer  content  of  the  beam,  which  changes  as  a  function  of 
total  flux,  does  not  affect  the  antimony  incorporation  kinetics. 
Ref.  [28]  has  shown  that  dimers  and  atoms  have  nearly  identical 
incorporation  probabilities  during  Group  V-rich  growth  and 
proposed  that  this  occurs  because  most  antimony  atoms  rapidly 
form  dimers  on  the  surface. 
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tetramers  [31,32],  and  ab  initio  cluster  model  re¬ 
sults  support  this  mechanism  [33].  In  this  mecha¬ 
nism,  two  atoms  from  each  tetramer  are  released  to 
the  gas  phase,  so  that  the  incorporation  coefficient 
is  never  larger  than  50%.  In  contrast,  100%  incor¬ 
poration  is  possible  with  As2.  These  stoichiomet¬ 
ries  have  been  widely  assumed  for  III-V  MBE 
growth  since  the  early  investigations  of  GaAs 
growth  using  modulated-beam  mass  spectrometry 
[31,34],  however,  other  studies  have  reported  As4 
incorporation  coefficients  greater  than  50% 
[35-37],  Studies  of  antimonides  report  that  Sb4 
incorporation  coefficients  are  closer  to  100% 
[38,39], 

The  vertical  offsets  between  the  dimer-  and  tet- 
ramer-phase  boundaries  in  Fig.  4  (s4/g4  vs.  s2/g2) 
reflect  differences  in  the  relative  efficiency6  of  the 
two  Group  V  species  in  stabilization  vs.  growth 
(s4/s2  vs.  t/4/t/2  )-7  For  example,  the  GaSb  phase 
boundaries  for  Sb4  and  Sb2  are  nearly  coincident 
(s4/s2  x  04/02),  suggesting  that  fluxes  of  dimers 
and  tetramers  have  roughly  the  same  relative  ef¬ 
ficiencies  in  both  stabilization  and  growth  of  GaSb. 
In  contrast,  the  large  offset  between  the  AlSb  phase 
boundaries  implies  that  proportionally  more  Sb4  is 
needed  for  stabilization  than  for  growth  (s4/s2  > 
g^/g2)-  Unless  the  different  behaviors  of  AlSb  and 
GaSb  are  solely  a  consequence  of  stabilization  dif¬ 
ferences,  then  this  result  implies  that  one  or  both  of 
these  materials  has  growth  kinetics  for  Sb2  and  Sb4 
that  differ  from  the  model  proposed  for  As2  and 
As4  on  GaAs  [31,32]. 


6  Here  we  define  efficiency  as  inversely  proportional  to  the 
flux  required  for  a  given  process.  The  stabilization  flux  sx  ranges 
from  zero,  where  no  Group  V  species  desorb  from  ihe  surface 
(and  thus  do  not  need  replacing),  to  infinity,  where  no  Group 
V  species  stick  to  the  surface.  The  corresponding  efficiency 
ranges  from  infinity  to  zero,  respectively.  The  values  gx 
[Torr/(ML/s)j  are  proportional  to  the  Group  V  flux  required  for 
V-limited  growth,  per  unit  of  growth  rate,  and  range  from  one 
(when  Groups  V  and  III  atomic  fluxes  are  equal)  to  infinity.  The 
analogous  efficiency  ranges  between  one  and  zero,  respectively. 

7  The  ratios  s4/s2  or  34/32  by  themselves  are  not  meaningful 
numbers,  because  they  contain  unknown  (and  different)  ioniz¬ 
ation  efficiencies  for  dimers  and  tetramers.  In  ihe  comparison 
s4/s2  •>-►34/32,  however,  the  ionization  efficiencies  cancel  out  on 
left  and  right,  so  that  we  may  compare  the  ratios  for  stabilization 
vs.  growth. 


In  general,  the  ratios  s4/s2  and  g4/g2  need  not  be 
equal,  because  they  relate  to  different  processes. 
The  value  gx  involves  Group  V  atoms  that  perma¬ 
nently  incorporate  into  a  growing  crystal.  There  are 
frequent  interactions  between  III  and  V  species, 
and  a  large  fraction  of  Group  V  species  are  buried 
by  growth  before  they  can  desorb.  The  surface  is 
highly  disordered  and  has  a  high  density  of  island 
edge  sites.  For  (001)  surfaces,  gx  increases  weakly 
with  temperature  under  typical  growth  conditions 
[30],  presumably  as  more  Group  V  species  desorb 
prior  to  incorporation.  On  the  other  hand,  sx  re¬ 
lates  to  the  replacement  of  Group  V  atoms  that 
desorb  from  sites  of  the  prevailing  surface  recon¬ 
struction  and  increases  strongly  with  temperature 
(Figs.  3  and  4).  The  surface  is  flatter,  more  ordered, 
and  has  fewer  types  of  adsorption  sites.  This  strong 
distinction  between  stabilization  and  growth  ap¬ 
plies  to  those  processes  for  AlSb  and  GaSb  that  are 
considered  here. 

For  InAs,  there  is  a  weaker  distinction  between 
sx  and  gx,  and  a  comparison  is  less  meaningful. 
Here,  both  arsenic  stabilization  and  arsenic-limited 
growth  occur  near  the  phase  boundary  between 
As-  and  In-rich  reconstructions  [(2  x  4)  and  (4  x  2)]. 
In  fact,  the  arsenic  flux  corresponding  to  [gx  x 
(indium  flux)]  and  the  flux  sx  plotted  in  Fig.  3 
represent  points  on  two  different  indium  flux  iso¬ 
bars  (with  finite  and  zero  indium  flux,  respectively) 
of  a  single  flux-temperature  phase  diagram  for  InAs 
growth.  Ref.  [23]  discusses  such  phase  diagrams  for 
InAs.  The  comparison  is  complicated  further  by  the 
fact  that  the  InAs  phase  boundaries  in  Figs.  3  and 
4  were  obtained  for  temperatures  considerably 
above  where  gx  was  measured  ( ~  420°C).  Measur¬ 
ing  both  growth  and  stabilization  around  420°C  is 
impractical,  because  extremely  small  arsenic  fluxes 
would  be  needed  for  stabilization,  while  at  higher 
temperatures  (>500°C)  the  isobars  merge,  and  the 
distinction  between  growth  and  stabilization  is  lost. 

The  slopes  of  the  phase  boundaries  reflect  kinetic 
differences  in  the  ability  of  dimers  and  tetramers  to 
stabilize  the  surface  at  the  phase  boundary.  For 
a  surface  in  equilibrium  with  vapor,  the  composi¬ 
tions  of  incident  and  desorbing  fluxes  would  be 
equivalent,  and  the  slope  of  a  log  pressure  vs.  in¬ 
verse  temperature  plot  is  proportional  to  the  free 
energy  of  desorption  or  adsorption.  The  situation 
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in  this  work  is  better  described  as  “steady  state”, 
because  the  incident  Group  V  composition  (dimer 
and  tetramer  concentrations)  that  is  imposed  by  the 
experiment  will  generally  not  match  the  desorbing 
composition.  In  this  case,  the  slopes  and  offsets  of 
the  reconstruction  phase  boundaries  may  depend 
on  the  energies  and  widths  of  various  kinetic  bar¬ 
riers.  For  In  As  and  AlSb,  the  two  lines  have  differ¬ 
ent  slopes,  implying  that  the  kinetics  of  dimer  and 
tetramer  incorporation  are  different.  The  parallel 
lines  for  GaSb  suggest  that  there  is  a  common 
rate-limiting  incorporation  step  for  dimers  and  tet- 
ramers.  One  possibility  is  that  GaSb  has  only 
a  small  barrier  for  breakup  of  surface  Sb4  into 
dimers,  so  that  dimer  and  tetramer  kinetics  are 
dominated  by  the  barrier  for  subsequent  dimer 
incorporation. 

The  adsorption  of  Group  V  species  onto  III-V 
surfaces  is  influenced  by  the  surface  structure  and 
stoichiometry.  Unlike  the  antimonides,  InAs  has 
a  Ill-rich  reconstruction  on  the  high  temperature 
side  of  the  transition,  so  it  is  not  surprising  that  the 
qualitative  features  of  its  phase  boundaries  are  dif¬ 
ferent  than  at  least  one  of  the  antimonides.  How¬ 
ever,  the  similar  dimer  and  tetramer  phase 
boundary  slopes  for  GaSb  stand  out  compared  to 
AlSb,  considering  that  the  two  materials  have  sim¬ 
ilar  surfaces.  Above  the  phase  transition,  both  ex¬ 
hibit  a  (4x3)  bilayer-Sb  reconstruction  with 
a  small  fraction  of  Group  III  atoms  ([40],  see 
footnote  3).  Below  the  phase  transition,  AlSb  has 
a  c(4  x  4)  reconstruction,  while  GaSb  shows  (2  x  5)- 
like  reconstructions.  The  (2  x  5)  surfaces  have  the 
highest  Group  V  content  of  any  III-V  zincblende 
material  and  are  conducting  rather  than  semicon¬ 
ducting  [41].  These  unusual  features  may  contrib¬ 
ute  to  the  difference  in  AlSb  and  GaSb  reactivity 
with  the  two  types  of  antimony. 


5.  Conclusion 

The  surface  reconstruction  phase  diagrams  re¬ 
ported  here  should  serve  as  a  guide  for  optimizing 
growth  conditions  and  as  a  means  for  calibrating 
substrate  temperature  for  MBE  growth  of  the  6.1  A 
compound  semiconductors.  The  phase  boundaries 


for  GaSb  under  dimer  and  tetramer  antimony  flux 
have  similar  slopes  and  little  offset,  in  contrast  to 
those  for  AlSb,  suggesting  fundamental  differences 
in  the  kinetics  of  antimony  stabilization  of  the  two 
material  surfaces. 


Appendix 

In  order  to  apply  the  transmission  thermometry 
technique  to  a  given  substrate  material,  a  quantitat¬ 
ive  calibration  curve  is  required.  This  curve  gives 
the  inflection  point  of  the  optical  transmission 
spectrum  as  a  function  of  the  substrate  temper¬ 
ature.  Because  the  routinely  used  calibration  curve 
[14]  for  GaAs  substrates  is  known  to  be  inaccurate 
by  as  much  as  30°C,  we  present  here  updated  calib¬ 
ration  curves  for  semi-insulating  and  n  + -doped 
GaAs. 

The  technique  used  to  obtain  the  calibration 
curves  is  described  in  detail  in  Ref.  [16].  Briefly, 
a  vacuum  furnace  was  used  to  heat  the  GaAs  sub¬ 
strate,  and  a  spring-loaded  thermocouple  was 
placed  in  direct  contact  with  the  wafer.  Transmission 
spectra  were  recorded  for  a  series  of  substrate  tem¬ 
peratures.  In  Fig.  5,  the  inflection  point  of  the  trans¬ 
mission  edge  is  plotted  as  a  function  of  substrate 


GaAs  Substrate  Temperature  (°C) 

Fig.  5.  Calibration  curve  for  transmission  thermometry  of 
semi-insulating  and  n  +  -doped  (n  =  1018 cm-3)  GaAs.  The 
dashed  curve  represents  the  previous  calibration  from  Ref.  [14] 
for  both  types  of  substrates. 
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temperature.  The  fits  are  fourth-order  polynomials 
with  the  form  T  =  a  +  bE  +  cE2  +  dE3  +  eE 4, 
with  T  in  °C  and  E  in  cm-1.  The  coefficients  are 
a  =  -  9070.5,  b  =  4.8124,  c  =  -  8.1385  x  10“4, 
d  =  5.8157  x  10-8,  e  =  —  1.5514  x  10-12  for  semi- 
insulating  GaAs  and  a  =  —  2.6712  x  104,  b  = 
11.759,  c  =-  1.8343  x  10“ 3,  d  =  1.2436  x  10"7, 
e  =  —  3.1536  x  10“ 12  for  n+  GaAs.  For  these  cal¬ 
ibrations,  wafer  thickness  was  470  +  15  pm 
(SI-GaAs)  and  440  ±15  pm  (n  +  -GaAs).  The 
n+  doping  level  was  (1-1.5)  x  1018cm-3. 
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